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5 ' Abstract 

I The time evolution of black holes involves both the canonical equations of quan- 

I ' turn gravity and the statistical mechanics of Hawking radiation, neither of which 



(N 



contains a time variable. In order to introduce the time, we apply the semiclassical 
approximation to the Hamiltonian constraint on the apparent horizon and show 
that, when the backreaction is included, it suggests the existence of a long-living 
\^ • remnant, similarly to what is obtained in the microcanonical picture for the Hawking 

! radiation. 
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A longstanding issue in theoretical physics is how to quantize Einstein grav- 
ity. Since in the theory there are constraints whose algebra contains structure 
functions (which depend on the space-time coordinates), one can conceive (in- 
5^ \ equivalent [1]) manners of lifting the Hamiltonian and momentum constraints 

to quantum equations, thus hindering the formal solution to the problem. On 
a more physical ground, one might notice that it is relevant to have a quantum 
theory of gravity at our disposal only for systems with very strong (Planck 
size) gravitational fields, such as those one expects in the early stages of the 
Universe. Because of Hawking's discovery of black hole evaporation [2], one 
also expect that quantum (or semiclassical) gravity plays a role in determining 
the dynamics of the late stages in the life of collapsed objects. 

In a general situation, one has to deal with the infinite number of degrees of 
freedom of the gravitational field configurations which are physically distin- 
guished only modulo coordinate transformations. This gives the constraints 
the form of functional differential equations, for which there is little hope of 
finding general solutions, and one then tries to reduce the number of degrees 
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of freedom from the onset. One way of implementing this program is the back- 
ground field method [3,4], in which one assumes a certain background metric 
manifold and does perturbation theory on it (although it turns out to be non- 
renormalizable). Another method is to impose a large (space-time) symmetry 
prior to quantization, resulting in a so called mini- (or midi-) superspace for a 
finite (or discrete) number of canonical variables [5,6]. In the latter approach 
to quantum gravity there is no time and one expects that a time variable can 
be identified by resorting to the semiclassical approximation [7]. 

Unfortunately, there are systems which do not allow for such a reduction, and 
(spherical) black holes fall among these cases, because one cannot find a finite 
set of variables in which to express both the Hamiltonian and the momentum 
constraints consistently over the whole space-time manifold [8]. Years ago, 
Tomimatsu [9] proposed an alternative scheme of reduction by considering 
the support of the constraints just near the apparent horizon and quantizing 
the residual degrees of freedom. 

Hawking radiation can also be studied from the point of view of the statis- 
tical mechanics by assuming the horizon area ^ is a measure of the internal 
(microscopical) degeneracy a of black holes according to [10] 



for Hawking particles to tunnel out through the horizon. When the space- 
time is asymptotically fiat, it is known that the canonical ensemble is formally 
inconsistent and one should instead use the microcanonical ensemble in order 
to ensure global energy conservation in the system [11]. For large black holes 
the two descriptions are practically equivalent, but they then disagree when 
the mass becomes small (of order the Planck mass). Within this approach 
the detailed form of the metric (including the backreaction) affects the grey- 
body factor [13] and the constant c in Eq. (2) cannot be determined without a 
complete knowledge of the (non-local) physics near the horizon. In this respect, 
the microcanonical picture is complementary to the dynamical approach of 



The aim of this letter is precisely to study the consequences which can be 
derived from confronting the constraints of canonical (quantum) gravity lo- 

^ This relation suggests that black holes are extended p-branes [11] and was also 
obtained in string theory [12]. 
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and thus determines the probability 
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Ref. [9]. 
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calized on the apparent horizon as in Ref. [9] with the global microcanonical 
picture of the Hawking effect, with a particular emphasis on the role played 
by the backreaction of the emitted radiation. 

Let us start from the dynamical approach. Tomimatsu [9] considered a con- 
formally coupled scalar field in a four-dimensional spherically symmetric 
space-time with action [14] 
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The four-dimensional metric has been written as 



ds"^ = gab dx" dx^ + ip^ dQ'^ 



(4) 



where x"' = [t, r), dQ"^ is the line element of a unit two-sphere, (p = ip{t, r) the 
radius of a two-sphere of area AiTip^, gab = gab{t, r) a two-dimensional metric 
and R its scalar curvature. Adopting the ADM decomposition of a generic 
two-dimensional metric [15] 



gab 
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(5) 



with N the lapse function and the (radial) shift function, one obtains 
that the momentum and the Hamiltonian H coincide Q near the apparent 
horizon (defined by gu = 0), Hr = H/ \pi. This just leaves the Wheeler-DeWitt 
equation (Hamiltonian constraint) for the wavefunction = 



1 
2 



4m2 



2H, 



^ = 



where m is the black hole mass (in units with c = h = G 
in gu = —l + 2m/r. The classical momenta are given by 



(6) 



1) as it appears 



n™ = ^ (m + ^) , H^ = 4m2 0, (7) 

and a dot denotes derivative with respect to the time t measured in the ref- 
erence frame of a (distant) static observer. To summarize, the system near 

^ In a suitable gauge with 7 = 2 and its conjugate momentum 11^ = (/?/4. Conse- 
quently 7 is not dynamical [9]. 
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the apparent horizon has two physical degrees of freedom: the Bondi mass m 
computed on the outer surface of the apparent horizon (of radius (p = 2 m) 
and the scalar field which can be used to reproduce the outgoing Hawking 
radiation as we shall show below (see also Ref . [9] ) . 

Eq. (6) can be solved exactly by separation of variables, and one finds out the 
set of solutions 



*(m,0) = *oexph k4>\ \ , (8) 




where \l/o is a constant and k a continuous (complex) wave number. For k 
real one has "plane waves" corresponding to a static black hole, and for k 
imaginary Hawking's behaviour [2] is recovered [9]. Despite the simple form of 
the derivation, the above solution taken at face value is not very illuminating 
about the temporal evolution of the black hole, since there is no time depen- 
dence in the function in Eq. (8), and one needs, at least in principle, to build 
a suitable superposition of modes in such a way that a time variable can be 
defined (in terms of m and (p). 

In practice, one resorts to the semiclassical approximation. This has been 
attempted already in Ref. [9], but here we want to employ the more systematic 
Born-Oppenheimer approach (see, for instance, Ref. [7,16]). First of all, we 
assume that the complete wavefunction factorizes into "gravitational" {'ip) 
and "matter" (x) parts, 

^i^{m)x{(t>]m) , (9) 
from which we get, upon substituting into Eq. (6), the coupled equations 



(10) 

^mX= [Hep - {H^)\ X , 
where we have introduced the rescaled functions 

=-0e+^/(nm)dm ^ ~ ^^g-i/(fl^)dm ^ ^^^^ 

and the scalar field Hamiltonian 



Til ^ 
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In the above {A{m)) = J d0x*(0; m) A m) for any operator A. It is 
noticeable that Eqs. (10) do not contain terms corresponding to gravitational 
fluctuations of the form which is obtained when the momentum 11^ enters 
quadratically [16,17]. This signals the fact that, on restricting the support of 
the constraints just on the apparent horizon, one is actually freezing most of 
the quantum fluctuations for the system [18]. 

We now assume that the variable m behaves (almost) classically. Correspond- 
ingly, the gravitational wave function is a wave-packet peaked on a classical 
trajectory m = m{t), thus it is well approximated by the WKB form 



gi Sci [m] 
V dmScl 

where the classical action Sd in the exponent is evaluated along m{t). The 
time variable t can then be naturally related to the mass m by 



^.2 

dt 



_d_ 
dm 



m 



_d_ 
dm 



(14) 



To leading order (in £pmp — K) the flrst of Eqs. (10) then gives the Hamilton- 
Jacobi equation 



m 



4m^ 



(15) 



where we used 3^3 d — H^, and the second of Eqs. (10) is unaffected by the 
above approximation. If ( 11^ ) = k"^ and constant, we recover the result of 
Ref . [9] . In particular, lor k — i K2 one has 



which is the well known Hawking's law of black hole decay [2] implying a finite 
time of evaporation (as measured by a static observer). The numerical coeffi- 
cient K,2 is arbitrary and its value could possibly be computed by considering 
the constraints over the complete space-time manifold. 

Instead of assuming that the state of the scalar field remains a (momentum) 
eigenstate, we can use the second of Eqs. (10) to determine the evolution of 
the state of 0, thus including the backreaction of the change in mass onto the 
radiation. Upon introducing the time according to Eq. (14) and making use 
of Eq. (15), we get the Schrodinger-like equation 
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^^ = -^n,x.-H,x., (17) 

where Xs = X<^~' ^'^'^'^ 

Schrodinger equations with time-dependent parameters can be solved by in- 
troducing invariant operators [19], whose form is exphcitly known for the gen- 
erahzed harmonic oscillator (see Ref. [20] for the details). In order to apply 
this method to Eq. (17) with an arbitrary m = m{t), one must first "regular- 
ize" the Hamiltonian Ti.^ by generating a non- vanishing "oscillator frequency" . 
This can be done by introducing a (small) harmonic-oscillator-like potential 0, 



n^^n^+'-<p\ (18) 

where e — > at the end of the computation. One then finds the spectrum 
{\ n,t) \ n E IN} of (exact) solutions to Eq. (17) and, on using Eq. (15), 



= 2 {n,t I \ fiji) 
1 
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where the auxiliary function x = x{t) must satisfy 



(19) 



X \ m 



A^(m'-]-ex= ^ ": „ , (20) 
dt \ m J 4 m^ x'^ 

with suitable initial conditions. The latter must be derived from some physical 
argument, as we shall attempt in the following. We further note that, for 
an evaporating black hole, m < and the right hand side of Eq. (19) is 
consistently positive. 

We shall assume that, for large mass (m ^ 1), Hawking's behaviour (16) be 
recovered and determine the corresponding x = xuit) by expanding in powers 
of 1/m, 



xh = E4> (21) 

where the coefficients Xi = Xi{t). Inserting the sum (21) into Eqs. (19) and 
(20) yields, in the limit for e — 0, 

^ The e-term can be viewed as a small mass for the scalar particle (IR cut-off) or 
prescription to discretize the spectrum. 
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Fig. 1. Luminosity of a radiating black hole in the canonical (dashed line) and 
microcanonical (continuous line) pictures. 

and one can further check that Eq. (16) is a consistent solution to the system 
of Eqs. (19) and (20) to even higher orders (details will be given elsewhere 
[21])- 

Prom the microcanonical description we know that the law (16) must be cor- 
rected for small masses. In fact the correct number density of Hawking quanta 
is better approximated by the microcanonical expression [11] 



^ ^ exp(47rm2) ^ ^ 

where [[X]] denotes the integer part of X. In the limit m oo, N equals the 
canonical ensemble number density (Planckian distribution). The decay rate 
(the opposite of the luminosity) for a black hole is then given by 



oo 

m = E J^k (u;) Tkioj) dw , (24) 

k 

where F is the grey-body factor and the sum is over particle species and angular 
momentum [22]. The comparisons between the luminosity and time evolution 

of m in the canonical and microcanonical pictures are displayed, respectively, 
in Fig. 1 and Fig. 2 for the simplified case F = A; = 1. The basic feature of the 
microcanonical evolution is the existence of a long tail due to the suppression 
of (relatively) high energy modes. 
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Fig. 2. Time evolution of the radiating black hole in the canonical (dashed line) and 
microcanonical (continuous line) pictures. 

Bearing in mind the above result, we can try and perturb Eq. (16) by assuming 
an expansion in powers of 1/m, 



Aw? i^^fn^ 

(25) 

where ki — Hi{t) and yi — yi{t). It turns then out that Eqs. (19) and (20) 
admit as solutions 



- = + - + (26) 

with a > 3, Ka constant, and a corresponding function x — x{t) (whose form 
is of no relevance here [21]). Again the value of the coefficients in the above 
equation cannot be fixed within the dynamical approach. 

One can fix the values of K2 and Ka in such a way that the luminosity given by 
Eq. (26) approximates the microcanonical luminosity of Fig. 1 and the history 
of the mass approximates that obtained from the microcanical picture, at least 
for m > 1. In Fig. 3 we plot m = m{t) for kq = 10~'^ k,2 and compare it with 
the microcanonical history. Since we used the large m approximation in order 
to solve Eqs. (19) and (20), the plot is not reliable for m < 1, although it 
is certainly suggestive that both curves show the presence of a long-living 
remnant of sub-Planckian mass. 



We wish to conclude by remarking that the main result of this letter is given 



8 



0.25 



1.75 



0.75 



1.25 



1.5 



0.5 




2000 



4000 



6000 



8000 



10000 



t 



Fig. 3. Time evolution of the radiating black hole from the Hamiltonian constraint 
(continuous line) as compared with the microcanonical picture (dashed line) for the 

case described in the text. 

by the set of coupled Eqs. (19) and (20) which govern the dynamics of an 
evaporating (apparent) horizon and allow one to study the backreaction of 
the emitted radiation. Unfortunately, they are rather complicated and do not 
allow for a straightforward analytical treatment. Further, one needs to spec- 
ify the initial conditions in such a way that relevant physical situations are 
described. This far, we have just investigated that system of equations by em- 
ploying the large mass approximation and compared with the results obtained 
from the statistical mechanics of the Hawking radiation in order to include the 
backreaction. We hope to carry on a more complete (possibly numerical) anal- 
ysis and gain a better understanding of the initial conditions in a forthcoming 
work [21]. 

References 

[1] P. Hajiceck, Nucl. Phys. B (Proc. Suppl.) (2000) 88 114. 

[2] S.W. Hawking, Nature 248 (1974) 30 and Comm. Math. Phys. 43 (1975) 199. 

[3] B. S. DeWitt, in Relativity, groups and topology, Vol. 1, DeWitt C. and DeWitt 

B. S. editors (Gordon and Breach, New York, 1964). 

[4] N. D. Birrell and P.C.W. Davis, Quantum fields in curved space (Cambridge, 
Cambridge University Press, 1982). 

[5] B. S. DeWitt, Phys. Rev. 160 (1967) 1113. 

[6] J. A. Wheeler, in Batelle rencontres: 1967 lectures in mathematics and physics, 

C. DeWitt and J. A. Wheeler editors (Benjamin, New York, 1968). 

[7] R. Brout and G. Venturi, Phys. Rev. D 39 (1989) 2436. 



9 



[8] CO. Lousto and F.D. Mazzitelli, Int. J. Mod. Phys. A 6 (1991) 1017. 
[9] A. Tomimatsu, Phys. Lett. B289 (1992) 283. 

[10] J.D. Bekenstein, Phys. Rev. D 7 (1973) 2333; Phys. Rev. D 9 (1974) 3292. 

[11] B. Harms and Y. Leblanc, Phys. Rev. D 46 (1992) 2334; Phys. Rev. D 47 (1993) 
2438; R. Casadio and B. Harms, Phys. Rev. D 58 (1998) 044014; Mod. Phys. 
Lett. A17 (1999) 1089. 

[12] A. Strominger and C. Vafa, Phys. Lett. B379 (1996) 99; J. Maldacena, Black 
holes in string theory, |hep-th/9607235 . 

[13] D. Page, Phys. Rev. D 13 (1976) 198; Phys. Rev. D 16 (1977) 2401. 

[14] P. Thomi, B. Isaak and P. Hajicek, Phys. Rev. D 30 (1984) 1168; P. Hajicek, 
Phys. Rev. D 30 (1984) 1178. 

[15] C.W. Misner, K.S. Thorne and J. A. Wheeler, Gravitation (San Francesco: 
Freeman, 1973). 

[16] C. Bertoni, F. Finelh and G. Venturi, Class. Quantum Grav. 13 (1996) 2375. 
[17] R. Casadio, Int. Jour. Mod. Phys. D 9 (2000) 511. 
[18] F. Englert, Phys. Lett. B 228 (1989) 111. 

[19] H.R. Lewis and W.B. Riesenfeld, J. Math. Phys. 10 (1969) 1458. 
[20] X.C. Gao, J.B. Xu and T.Z. Quian, Phys. Rev. A 44 (1991) 7016. 
[21] R. Casadio, in preparation. 

[22] F. Halzen, E. Zas, J.H. MacGibbon and T.C. Weeks, Nature 353 (1991) 807. 



10 



